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ABSTRACT 

A common characteristic of the interaction between the coastal topography and 
eastern boundary currents (EBC) is the appearance of cold filaments and mesoscale eddies. 
Hydrographic and satellite temperature data obtained during a cruise on board R/V Point 
Sur off Point Arena, California, in May 1993 were analyzed to study a particular eddy field 
in this area. 

The hydrographic data was first used to verify the remotely sensed surface temperature 
field, using three dimensional data visualization. Selected vertical levels from each 
hydrographic station were then interpolated into a broader, finer resolution grid domain in 
preparation for an eventual model initialization, using multiquadric interpolation. 

The results verify the existence of an eddy and show its signature in the vertical to 
about 300 meters depth. A sensitivity study of interpolation parameters was performed to 
evaluate approximately the optimal set of parameters, showing that the multiquadric 
interpolation resolves very well the temperature field in the upper levels and introduces 
small amplitude, small scale noise in the deeper levels. This noise can be eliminated by a 
more thorough parameter sensitivity study. 

A comparative analysis of satellite AVHRR data, in situ hydrographic data and 
interpolation results shows that the eddy did not experience significant horizontal 
translation during the period of the cruise suggesting that the asynopticity between the 
different data sets is not significant for the purpose of objectively analyzing the temperature 
field. 
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I. INTRODUCTION 



Following the early interest in the strong and well developed western boundary 
currents such as the Gulf Stream and the Kuroshio, the oceanographic community is 
becoming aware of the importance of the eastern boundary currents in the dynamics of the 
ocean. This is particularly true in the coastal zone, where the interaction between the 
eastern boundary currents and the boundary is strongest. 

Eastern boundary currents (EBC) are typically characterized as a slow, broad and 
shallow equatorward flow with diffuse boundaries or transition zones separating them from 
the local coastal currents. Gill (1982) summarizes several aspects of EBC, including 
driving mechanisms, non-local influences due to coastally trapped waves, bottom friction 
and topography effects, upwelling, frontogenesis and turbulence. An intricate circulation 
results from this complex system where cold filaments extending offshore and permanent 
and recurring eddies are common features. 

Eddy activity in EBC regions has been observed as early as the mid nineteen sixties, 
but the general perspective of the extent of eddy phenomena in time and space was only 
made possible by the more recent availability of satellite observations (Pares-Sierra et al. y 
1993). Since then, eddies and filaments have been observed in eastern ocean boundaries 
using different types of remotely sensed and/or in situ observations (e.g., Pares-Sierra etal., 
1993; Sousa, 1992; Strub etal., 1991; Ramp etal., 1991; Haynes etal., 1991). 

A useful technique for the study of eddies and cold filaments is numerical simulation 
(e.g. Pares-Sierra etal . , 1993; Batteen et al., 1992). Numerical models require some form 
of observed or statistical data for initialization, forcing and verification. Although satellites 
are becoming an important source of regular and nearly synoptic observations, 
oceanographic data is in general sparse, and spatial and temporal gaps are common, 
presenting difficulties in model initialization or data assimilation. 
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Several methods to interpolate and optimize observed fields for use in models have 
been developed (see Haltiner and Williams, 1980, for a good description of the interpolation 
and initialization methods commonly used in numerical modeling). These methods use data 
correlation and/or historical data to relate observations collected in an irregular 
observational space to information projected into a regular model grid. The transformation 
must be able to keep the interpolated fields within realistic bounds and be able to preserve 
the state of slight imbalance between the mass and motion fields that ultimately drives the 
system. The purpose of this study is to objectively grid and analyze a coastal ocean eddy 
seen in satellite Advanced Very High Resolution Radiometer (AVHRR) and in situ 
hydrographic data, making use of the fairly recent multiquadric interpolation method 
(Hardy, 1971; Nuss, 1993) to interpolate and prepare observed fields in two and three 
dimensions. Multiquadric interpolation is the chosen method because it has proven superior 
qualities over other methods (e.g., Barnes and Cressman objective analyses) in the 
objective analysis of the observed fields and reduction of spurious high wavenumber noise 
in space (Nuss and Titley, 1993). 

Data collected from 21 to 25 May 1993 by R/V Point Sur off Point Arena, California, 
is processed with the ultimate objective of being used to initialize a coastal ocean 
circulation model. We begin by describing the data and its collection methods in the next 
section, followed by data analysis, multiquadric interpolation description, and results. The 
final section includes conclusions and suggestions for future work. 
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II. DATA COLLECTION 



A. SATELLITE DATA 

The satellite data used in this work consists of images obtained from the Oregon State 
University. The images were from channel four of the Advanced Very High Resolution 
Radiometer (AVHRR) flying on the National Oceanic and Atmospheric Administration 
(NOAA) series of Earth observing satellites. This sensor observes the surface of the Earth 
in the infrared region (IR) with wavelengths between 10.5 and 11.5 |im. The temperature 
data were computed by converting radiometer counts to temperatures directly without any 
correction for water vapor (atmospheric correction), because no other channels were 
available to perform a multispectral correction. 

After downloading the images from Oregon State University via Internet, further 
processing was done at the Naval Postgraduate School using software developed at the 
University of Florida, Miami. This processing consisted of correcting the images for 
appropriate navigation (geographical correction) and remapping the entire image into the 
region of interest, centered at a point off Point Arena, California. The ground resolution 
depends on the distance between each particular pixel and the satellite nadir and is in the 
range of 1 to 3 km. A suitable colormap was then applied to enhance the range of 
temperatures of the ocean surface and make more visible the regions contaminated by 
clouds. A mask was applied over land. 

The satellite derived temperatures were obtained by forcing a linear relationship 
between pixel value and temperature value. To obtain the satellite temperature series (SAT) 
in Figure 2, for example, the image for day 22 May 1993 was sampled at the same positions 
of the CTD stations, and the pixel value was read and divided by ten to convert to 
temperature, thus providing a resolution of one tenth of a degree Celsius. 
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Most of the images collected during the period of the cruise contained large amounts 
of cloud cover; for this reason, only the images for days 21 and 22 May (23 Greenwich 
Mean Time) were used in this work (see Figures 3 and 4). 

B. HYDROGRAPHIC DATA 

The hydrographic data used in this study was collected during a Naval Postgraduate 
School cruise on board R/V Point Sur, between the 21* and the 25* of May, 1993, using a 
Neil Brown Mark III Conductivity, Temperature and Depth (CTD) instrument. The CTD 
casts were performed on an array of stations listed in Table 1 (see also Figure 1). 



TABLE 1; CHARACTERISTICS OF THE CTD STATIONS 



LONG (W) 


LAT(N) 


NUMBER 


DATE 


TIME ( h ) 


MAXZ(m) 


-123.697 


38.546 


1 


22 


23.9 


500 


-123.779 


38.643 


2 


22 


22.7 


460 


-123.857 


38.750 


3 


22 


21.4 


410 


-123.951 


38.859 


4 


22 


20.0 


455 


-124.027 


38.970 


5 


22 


18.5 


530 


-124.176 


38.889 


6 


23 


10.6 


2420 


-124.080 


38.782 


7 


23 


8.3 


2110 


-123.984 


38.675 


8 


23 


6.0 


1965 


-123.917 


38.587 


9 


23 


4.0 


1940 


-123.811 


38.478 


10 


23 


1.8 


1630 


-123.940 


38.393 


11 


23 


24.0 


2280 


-124.026 


38.493 


12 


23 


21.2 


3125 


-124.100 


38.591 


13 


23 


18.8 


2500 


-124.212 


38.700 


14 


23 


16.1 


3200 


-124.313 


38.808 


15 


23 


12.8 


3260 


-124.456 


38.703 


16 


24 


11.4 


3500 
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TABLE 1: CHARACTERISTICS OF THE CTD STATIONS 



LONG (W) 


LAT(N) 


NUMBER 


DATE 


TIME (h) 


MAXZ(m) 


-124.332 


38.595 


17 


24 


9.1 


3445 


-124.239 


38.507 


18 


24 


6.4 


3430 


-124.131 


38.414 


19 


24 


4.2 


3430 


-124.027 


38.318 


20 


24 


2.0 


3430 


-124.149 


38.236 


21 


25 


11.6 


3600 


-124.248 


38.338 


22 


25 


9.4 


3670 


-124.356 


38.422 


23 


25 


0.6 


3650 


-124.467 


38.525 


24 


24 


22.0 


3605 


-124.574 


38.621 


25 


24 


19.6 


3660 


-124.787 


38.807 


26 


24 


14.5 


3610 


-124.681 


38.714 


27 


24 


16.9 


3650 


-124.498 


38.356 


28 


25 


3.8 


3760 


-124.629 


38.282 


29 


25 


6.0 


3840 


-124.049 


38.150 


30 


25 


13.6 


3610 


-123.948 


38.061 


31 


25 


15.7 


3510 


-123.842 


37.966 


32 


25 


17.9 


3480 



The instrument’s pressure, temperature and conductivity sensors were calibrated prior 
and after the cruise. A General Oceanics rosette sampler was used to collect in situ water 
samples for conductivity calibration after the cruise. The CTD sampling rate was 32 Hz. 
The raw data was collected and processed using a software package developed at EG&G 
Marine Instruments. The processing done at this level consisted of removing bad data 
points (oudiers) using a first difference criteria to detect suspicious data points and linear 
interpolation to perform the correction. The data was then averaged to every 2 dbar and 
recorded in files for posterior data analysis. 
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An underway data acquisition loop was active during the cruise, automatically 
recording averages of temperature and salinity data at the 2 m depth (among other data) 
every 30 seconds (Ship’s Data Acquisition System). This surface data was used in this 
study to fill data void areas and will be hereafter referred to as SAIL data (see section 
m.B.2). 

Moored current meters (Acoustic Doppler Current Profiler - ADCP - and 
conventional rotor meters) were also deployed at station numbers 3, 7, 8, 9, 13, 17, 18, 19 
and 23. The ADCP is able to measure the three components of the current using an array of 
four acoustic beams; the Doppler shift induced in the returning echoes by the moving 
column of water is processed to compute the absolute three dimensional current vector. An 
appropriate time gating of the echo yields a profile of the velocity field in the vertical 
between the instrument and the bottom for downward looking ADCP sensors, or between 
the instrument and the surface for upward looking sensors, at several depth bins (for more 
details on the ADCP instrument, see Acoustic Doppler Current Profilers, 1989). 

The ADCP sensors deployed in this cruise were looking upward from the top of the 
mooring, providing current measurements from 300 m depth to the surface. In addition to 
the moored ADCP instruments, a hull mounted ADCP looking downward recorded 
velocity data continuously following the track of the ship. An example of the final product 
combining the moored and hull mounted ADCP data is presented in Figure 5. This figure 
shows the horizontal current vectors corresponding to the depth bin 15-23 m with velocities 
in cm s* 1 . This current meter data was used in the analysis of the surface temperature field. 
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in. DATA ANALYSIS 



A. SATELLITE DATA 

Since the satellite data consisted only of radiances measured by channel four of the 
AVHRR, the analysis of the satellite data collected during the cruise period presents two 
major problems: the comparison of the remotely sensed sea surface temperatures (SST) 
with in situ observations in order to determine the amount of atmospheric influence on the 
sensed image and the correction of the satellite data for atmospheric influences in order to 
make the data usable in later stages of the work. 

Figure 2 illustrates the method used to determine the corrections to be applied. The 
CTD series was obtained by averaging the first five meters of temperature readings at each 
station and ordering the results by time of CTD cast. The satellite series was obtained by 
sampling the image in Figure 4 at the positions of the CTD stations. These two spatial series 
contain a time lag introduced by the lack of synopticity of the CTD data; the time of the 
satellite image is comparable to the time of the first few CTD stations. The difference 
between the two series is thus not only due to atmospheric effects (which are not time 
dependent) but also due to a time-varying term (see Figure 2, middle plot). This latter term 
contains information on the movement of the eddy shown in Figure 4 that can be estimated 
by carefully comparing the CTD data collected in the final pan of the cruise with the 
satellite data from two or three days earlier. This information on the motion of the eddy 
during the period of the hydrographic cruise can be useful for adjusting the analysis for 
asynopticity, and to evaluate the performance of a prediction model started from the 
analysis. 

The problem of making the satellite observations usable can then be regarded in two 
ways: either making the CTD stations “more synoptic” relative to the satellite data or 
adjusting the satellite data to the weakly synoptic CTD data that we know is the “truth” in 
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its time/space frame, correcting the atmospheric term in the process. Both approaches 
remove the trend in the temperature difference in the middle frame of Figure 2, in a least 
squares sense, and shift the satellite series over the in situ series. This adjustment was 
accomplished by fitting least squares straight lines to both series, computing the 
point-by-point difference between these two lines and finally adding the point-by-point 
difference to the original satellite series. The bottom plot in Figure 2 shows that the results 
obtained are very satisfactory. A comparison of the variability of the corrected satellite 
series relative to the CTD series and the variability of the CTD series itself is given by the 
ratio: 



([ t ctd~ t sat1 ) 

{[T cw -{T ctd )} 2 ) 



(III.1) 



where Tcro is the CTD temperature series, and T w is the series of adjusted satellite-derived 
temperatures taken at the positions of the CTD stations. The angled brackets represent 
ensemble averages. For the data shown in the bottom frame of Figure 2 this ratio is 0.44. 
The mean square difference between the adjusted satellite series and the CTD series is 
thus smaller than the variance of the CTD series, meaning that the values of the adjusted 
satellite series may be regarded as representative of the sea surface temperature. 

This procedure guaranties the best possible match between the two series at each CTD 
station position and also inside a domain defined by the outer stations through bilinear 
interpolation of the correction terms. For positions outside this domain, the satellite 
readings may be corrected using averages of the nearest-neighbor correction terms. 

The purpose of these corrections is to be able to use satellite-derived sea surface 
temperatures with a minimum atmospheric effect in the multiquadric interpolation outside 
the area defined by the CTD stations, since it is desired to analyze a region somewhat larger 
than that defined by the CTD array (see Figures 4 and 9). The satellite-derived temperatures 
will influence the surface values for the areas not covered by CTD stations in the horizontal, 
and this will augment the data for interpolating in the vertical for the underlying column of 
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water. It will be shown later that this influence is significant in the final surface temperature 
maps produced by the multiquadric interpolation. 

B. HYDROGRAPHIC DATA 

1. Analysis of CTD Stations 

The data collected in the CTD casts as described in the previous chapter were first 
analyzed to determine to what extent the features visible in Figure 4 were not an 
atmospheric phenomenon, such as local wind swirls induced by coastal topography. For 
this purpose a visualization software package available from Silicon Graphics, Inc. was 
used (Explorer, 1993). 

The CTD stations provide both temperature and salinity data over a range of 
depths. For simplification in the visualization, only the stations lying within a more or less 
regular square were used. This includes stations 1 through 25. The maximum depth 
considered was 2001 m; station numbers 1 through 5 and 8 through 10 do not reach this 
depth due to the shoaling bottom. Although bottom topography will be included in future 
work, a similarity method was used here to extend the profiles below the shallow pans of 
the bottom topography down to the 2001 m depth, in order to achieve a regular three 
dimensional domain. This method searches for temperatures in the nearest profile equal to 
the deepest temperature in the profile to be extended. After the match is found, the lower 
part of the nearest profile is copied to the profile to be extended. The similarity method is 
attributed to Helland-Hansen, but no citation was found. The method reflects the theory that 
a horizontal temperature (or any other conservative property) line in a particular region can 
be found in the vertical in some position away from this region due to the general oceanic 
circulation (slow sinking and spreading). 

The vertical resolution of the raw CTD data is one data point every 2 m, starting 
at z=-l m. To reduce computer memory requirements, the profiles were sampled at 52 
preselected levels (every 2 m from -1 to -1 1, every 4 m from -11 to -19, every 6 m from - 
19 to -43, 8 m to -51, every 10 m from -51 to -101, every 20 m from -101 to -301, every 
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30 m from -301 to -481, 40 m to -521, every 50 m from -521 to -671, 60 m to -731, 70 m 
to -801 and every 100 m from -801 to -2001). 

The data were visualized by Explorer using a perimeter coordinate system, with 
uniform grid spacing in the horizontal (approximately 14 km) and arbitrary spacing in the 
vertical, to allow selection of different levels and the setting of an appropriate aspect ratio 
for visualization. Some images resulting from the CTD visualization are shown in Figures 
6 through 8. Of particular interest is Figure 8 showing a map of the surface temperatures as 
measured by the CTD. Since only 25 CTD stations were used in the visualization, the 
locations of the grid points are clearly evident due to the coarse resolution. However, a 
comparison of this map with Figure 4 shows that the main concentrations of cool and warm 
water were also defined by the CTD casts, thus verifying that Figure 4 does indeed 
represent an oceanic eddy. At the same time, Figure 7 shows a fluctuation of the main 
thermocline that corresponds well to the warm and cold waters at the surface, 
demonstrating the vertical manifestation of the eddy. The eddy appears to have a definite 
signature in the upper 300 m, in agreement with other sources of data including current 
meter data. Further discussion on the correspondence between the raw CTD data, the 
satellite images and the results of the interpolation scheme will be presented in Chapter IV. 

The vertical structure of the data is in fact quite complex. Isothermal surfaces are 
often found folded several times, particularly in the upper levels, perhaps due to 
observational errors. However, it is expected that the vertical motions in a coastal eddy like 
the one considered here may be large, and these could lead to vertical distributions of 
temperature and salinity that would allow folding of isotherms and still maintain the 
hydrostatic stability of the water column. Static stability analysis was not carried out in the 
present raw data analysis, but it will be considered in future work. Both temperature and 
salinity data are available for the computation of the three dimensional density field that 
can be visualized in the same way the temperature field is, providing a quick way to 
determine the stability of the entire water column defined by the CTD stations. 
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In order to improve the visual appearance of the raw CTD data, and to prepare it 
for input into an ocean prediction model with approximately 2 km resolution, the CTD 
observations must be interpolated onto such a grid. The method of multiquadric 
interpolation used for this purpose is described in the next section. 

2. Multiquadric Interpolation 

We begin with a description of the method, following closely that provided by 
Nuss (1993). Consider a three dimensional irregular surface containing a field sampled at 
N observation points. For each observation point we know the location (x,y,z) and the value 
of the field (//). This irregular surface can represent almost any scalar field; for this work, 
H is temperature. 

The interpolation of the field value in any arbitrary point X is based on a weighted 
sum of radial basis functions: 

N 

H(X) = £a,e(X-X,). (Ill .2) 

i = 1 

where cx, is the weight of the kernel function Q , its argument being the radial distance 
between the point X and the observation point X,. For multiquadric interpolation, the 
kernel function is a hyperboloid function: 

<2(X-X,.) = [ (X-X,) 2 + c 2 ] 1/2 . (III.3) 

The constant c is an arbitrary small number to keep the function from vanishing 
at the observation points. For three dimensional interpolation, the kernel function is 
therefore: 

<2(X-X ( .) = [ (x-x,) 2 + (y-y i ) 2 + (z-z,) 2 + c 2 ] 1/2 . (III.4) 

To solve for the weights, we apply Equation III.2 to all N observation points, 
resulting in a set of N equations with N unknowns: 
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(111*5) 



N 



H ( X j) = '£<* i Qi( x j>yj> z j) 



i = 1 



In matrix notation, the solution is given by: 




(III.6) 



where Q* is a square matrix of equation coefficients. This matrix was inverted using 
Matlab’s matrix inversion function (Madab, 1992), in the two dimensional case, and using 
routines available in the LINPACK library of mathematical routines (Unpack, 1978), in 
the three dimensional case. 

Once the weights are known, the interpolated solution to any point (e.g., regularly 
spaced grid points) is then given in matrix notation by: 



(x t ,y v z t ) being the coordinates of each regular grid point. The two dimensional problem is 
just a particular case of the more general three dimensional case, with all z coordinates set 
to zero. The matrix on the left hand side has as many rows as there are grid points, but the 
number of columns is always equal to the number of observation points. In other words, 
the solution depends only on the number of observation points available and is 
independent of the size of the grid to be interpolated. 

The method presented above is valid for noise-free observations. Unfortunately, 
real data include observational errors that must be taken into account. This is achieved 
through the introduction of an uncertainty parameter and a smoothing parameter (Nuss and 
Titley, 1993) in Equation III.5. In matrix notation: 



H 




(III.7) 



where: 




(III.8) 



H. = [Q ij + NXaf-5 iy ] a.. 



(HI*9) 
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where X is the smoothing parameter, or/ is a vector of uncertainties in the observations 
(mean-squared observation error) and 8,, is the Kronecker delta function. The difference 
now is that the diagonal elements of Q<, are modified to include a factor involving those 
two parameters. Note that the uncertainties vector can have different values for different 
observation points, an useful property that gives the user of the method enough freedom to 
optimally set the parameters for best fitting any part of the grid. 

There are several judgements that must be made when this method is to be used. 
The first two are selecting the values of the smoothing parameter X and the interpolation 
constant c. The smoothing parameter determines how damped the interpolated field will be, 
and its value will depend on the application and kind of data being interpolated; large 
values imply strongly damped analyses. The interpolation constant has also some 
smoothing associated with it, but its major effect is on the shape of the basis functions. 
Large values imply rather flat interpolating functions and thus tight gradients or closely 
spaced observation points are not well represented. The values for these parameters depend 
to a large extent on the particular application. While research is still going on to determine 
optimum values for c and X, setting them remains a trial and error process. 

Another factor to be considered is the aspect ratio of the vertical to the horizontal 
scales. In the ocean, as well as in the atmosphere, observations are much more closely 
spaced in the vertical than in the horizontal. This is needed because the scales of ocean 
variability are much shorter in the vertical than in the horizontal. Thus, the vertical 
coordinates must be expanded in such a way that the distances ( X-X t ) in Equation 113.4 are 
effectively nondimensionalized by horizontal and vertical length scales appropriate for the 
data being analyzed. The proper amount of vertical coordinates expansion must be 
determined and here again we are forced to use a trial and error process, although the aspect 
ratios in common practice in oceanography provide a guideline. The best results were found 
using an aspect ratio of approximately 1/2100 corresponding to normalizing the horizontal 
coordinates (dividing by the length of the shortest grid domain axis) and leaving the vertical 
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coordinates unchanged. Scaling the vertical dimension by dividing by the maximum depth 
(aspect ratio of approximately 1/44) yielded very poor results, but there should be several 
possible combinations of the smoothing parameter, interpolation constant, mean-square 
errors and aspect ratio that can give reasonably good results. For a more complete 
discussion of the sensitivity of the method, and procedures to help specify the scaling and 
smoothing parameters, the reader is referred to Nuss and Titley (1993). 

The multiquadric interpolation method outlined above was applied to the CTD 
stations after the processing described in the previous section and after appending four new 
“stations” with SAIL data to the vector of observations (see Figure 9). SAIL points were 
selected with a spacing similar to the CTD stations spacing, and they contain information 
on temperature and salinity only at the surface. The main purpose of including these four 
SAIL stations was to determine to what extent the lack of observations directly below 
surface observations influences the interpolated fields at depth (in this way simulating the 
future inclusion of satellite data in the interpolation). A minor draw back is that the 
inclusion of the SAIL data points increases the asynopticity problem, since they were 
measured on the 21" of May, whereas the CTD casts were collected between the 22^ and 
the 25 th of May. 

The multiquadric interpolation was applied in both two and three dimensions 
separately, with the purpose of comparing the surface temperature analysis produced by the 
three dimensional version with the analysis using the two dimensional version. The latter 
uses only the surface temperature values from CTD and SAIL on a 48x45 grid; the grid 
spacing is 1825 m in the x direction and 1831 m in the y direction. A rotation in the grid 
domain was applied so the x axis is approximately parallel to the station line with the 
greater number of stations running roughly in the NW-SE direction (see Figure 9). The 
coordinate system used is the Universal Transverse Mercator (UTM) with origin of 
coordinates in the lower left comer of the grid domain (see the U. S. Geological Survey 
Professional Paper #1395, 1987, for transformation formulae). The limits of the domain 
were arbitrarily set in order to include the eddy and the cold water in the upper left comer 
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but to exclude the warm filament to the south (see Figure 4). Some CTD and SAIL stations 
were intentionally excluded from the domain, also. These data points were included in the 
computation of the weights and helped fill in the areas without data. 

The three dimensional version interpolates grid points in the same horizontal 
layout as in the two dimensional application, but extending in the vertical down to -2001m. 
The same preselected levels introduced in section III.B.l were used here for observation 
points. The domain was then organized in a long vector of observations to compute the 
weights and perform the interpolation. The interpolated levels coincide with the 
observation levels, but this is not necessary in general. Several runs using different 
interpolation parameters and mean-square errors for the observations were performed, to 
make a rough estimate of the optimal parameters. The results of the two and three 
dimensional versions were then visualized using the Explorer visualization package. 

The results obtained with the multiquadric interpolation of the CTD and SAIL 
data are presented and discussed in the next chapter. 
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IV. RESULTS 



This chapter presents the results obtained with multiquadric interpolation of the raw 
CTD and SAIL data into the grid domain starting with the two dimensional case, which is 
then compared to the surface map generated from the three dimensional multiquadric 
results. 

A. TWO DIMENSIONAL MULTIQUADRIC INTERPOLATION 

The two dimensional multiquadric interpolation results are presented in Figure 10. 
These results were obtained by applying a two dimensional version of Equation III.9 to the 
four SAIL data points and to the average of the first five meters of each CTD station (first 
three levels), interpolating onto the grid described in the previous section. We are able to 
identify the major features of the temperature distribution at the surface seen in the two 
dimensional results by comparing with Figure 4. These features include the regions of 
warm water on the right side of the domain and at the lower left comer, as well as the cooler 
water in the upper part of the image. The method was also able to reproduce the general 
regions of cold and warm water in the interior of the domain that can be seen in the satellite 
image. Of course, an exact match of the interpolated results with the satellite image is not 
possible to achieve, in the first place because the satellite image is a high resolution 
snapshot of the surface temperature field at an instant in time near the beginning of the 
cruise, and in the second place because of the temporal bias between the remotely observed 
temperatures and the CTD observations, discussed in section III.A. 

The influence of the lack of synopticity between the satellite derived temperatures and 
the CTD and SAIL observations on the analysis can now be further quantified, with the aid 
of Figures 3, 4, 5 and 10. The low resolution of the CTD and SAIL sampling introduces a 
high degree of uncertainty in the identification and correspondence between these figures 
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of mesoscale features, such as the center of the eddy or other features in the outer limits of 
the eddy. 

The determination of the exact location of the center of the eddy on each of the Figures 
3, 4 and 10 gives an idea of the translation of the eddy in the horizontal, since these figures 
present data collected at different instants in time. The interpretation of features other than 
estimation of the location of the center of the eddy is difficult because there are several 
mechanisms that can contribute to a temporal change in the shape of the temperature field 
associated with the eddy (e.g., secondary circulation, horizontal advection, spinning). 
Furthermore, a positive identification of the existence of an eddy is somewhat uncertain. 
However, the location of the center of the eddy can be estimated with some precision from 
the satellite images in Figures 3 and 4 and the multiquadric interpolation results in Figure 
10. An analysis of Figure 5 shows that the center was located approximately at 124.3° W 
and 38.7° N during the one- week period of data collection. This result is supported by the 
satellite images, in particular Figure 4 that shows a spot of warm water in that same 
position, just to the northwest of station number 18. In the UTM grid, this station is 
positioned on grid lines x=34 and y=17, so the pocket of warm water northwest of the 
station is positioned approximately at grid points x=29 andy=16. Figure 10 shows a warm 
spot at this position and it is reasonable to assume that it represents the center of the eddy. 
The existence of a warm region at x=42, y=22 is the result of the strong influence of the 
relatively warm SAIL stations on the interpolation scheme, causing the generation of a 
local maximum in this area. 

From these combined analysis and observations, it appears that the eddy did not 
experience significant horizontal translation during the duration of the cruise. Therefore, 
although there is a time lag between satellite data and hydrographic data as shown in Figure 
2, it appears that the effects of time delay on the interpolation procedure are not significant 
because the eddy remained in approximately the same position. Figures 3 and 4 also show 
the eddy in the same position over a period of 24 hours, supporting this conclusion. Had the 
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eddy moved either alongshore or offshore, the discrepancies between the different data sets 
would be larger, and the level of uncertainty would be increased considerably. 

Nothing has been deduced about the rotational aspect of the eddy dynamics, because 
this is more difficult to detect from the available data. There is no evidence that the warm 
filament visible in Figure 4 that is located near station number 24 was sampled by this CTD 
cast (cf. Figure 8). This fact may have several interpretations. One might conclude that the 
feature moved due to rotation of the eddy. However, this may also be explained by 
secondary circulations or temporal changes in the field. 

The range of values in the grid domain is exactly the same as the surface temperature 
range observed in the CTD casts. However, the interpolated maximum and minimum 
temperatures do not occur at the same locations as the observed maximum and minimum 
temperatures because multiquadric interpolation, by definition, does not fit exactly with the 
observation points due to a non-zero interpolation constant and numerical truncations. 

With these considerations in mind, the method succeeds in keeping the interpolated 
values within realistic bounds and reproducing in fine detail the surface temperature 
structure. 

B. THREE DIMENSIONAL MULTIQUADRIC INTERPOLATION 

A more complete three dimensional multiquadric interpolation algorithm was used to 
interpolate the CTD observations onto the same horizontal grid domain and along the 
vertical down to -2001 m. Where necessary, CTD stations were extended as described in 
section m.B.l. Surface SAIL data points were also included in the input to the 
interpolation. The same levels mentioned in section III.B. 1 were used for vertical resolution 
of the interpolated field. 

A series of three dimensional interpolations were performed in a trial and error process 
in order to roughly estimate the best set of interpolation parameters, the smoothing 
parameter the interpolation constant c and the uncertainty level a, from Equation III.9. 
The initial values of these parameters were taken after Nuss (1993). For the purpose of a 
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first level sensitivity study, several combinations of parameter values were used. Table 2 
lists the values used in this work. 



TABLE 2: INTERPOLATION CONSTANTS 



CONSTANT 


VALUES USED 


X 


0.025 and 0.01 


c 


0.005,0.05 and 0.01 


a 


0.1,0.25 and 0.5 



A total of six combinations were used, comparing results with the observed CTD 
temperature values at the first three levels, at the last level and at the level of maximum 
temperature (this level is z=-3 m for CTD observations, and varied between z=-3 m and 
z=-15 m for interpolated values). The best set of parameters thus obtained was ^=0.025, 
a=0.25 and c=0.005. The three dimensional multiquadric interpolation results for this set 
of values are presented in Figures 11 through 13. 

Figure 1 1 shows the surface temperature map. Notice the close similarity between this 
surface map and the one in Figure 10, obtained using only the observed temperature values 
at the surface. The surface map obtained using the three dimensional algorithm does a very 
good job in maintaining the interpolated values in the same range and keeping the positions 
of the main features of warm and cold water, including the assumed center of the eddy. 
Figure 12 shows several isothermal contours and the isothermal surface of 8 °C. The 
primary purpose of this figure is to illustrate two results of the three dimensional 
interpolation: the resulting smooth isothermal surface and the close agreement of this 
surface with the same surface produced by the analysis of raw CTD data (Figure 6), 
showing that with the chosen set of parameters the method produces good results both in 
the horizontal and in the vertical. 

The three dimensional results, as far as smoothing is concerned, can be affected not 
only by the smoothing parameter X and uncertainty level a, but also by the vertical 
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distribution of observation points. As mentioned previously, the vertical levels of observed 
points were selected to give more resolution in the upper ocean and less resolution as depth 
increases. It was noticed that this fact is probably the cause of the appearance of small scale 
noise in the interpolated deeper levels, mostly below z=-1000 m. This is evident in Figure 
13. The main feature of this figure is the isothermal surface of 3 # C with small amplitude, 
small scale ripples. These ripples are artificially introduced by the interpolation algorithm. 
A possible way to reduce this noise is by setting the <J parameter as a vertically varying 
parameter, keeping the level of smoothing in the upper levels and increasing the smoothing 
in the lower levels. Another possibility to reduce the noise is to make a thorough parameter 
sensitivity study, using reference criteria other than the one mentioned above, or extending 
that criterion to include the resulting variability in the deeper levels. This way, a different 
optimal set of parameters could be selected. Both approaches will be taken in future work. 
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V. CONCLUSIONS AND RECOMENDATIONS 



Hydrographic and satellite-derived data were analyzed and processed to evaluate them 
for possible input in a coastal circulation model. The hydrographic data consisted of a set 
of 32 CTD casts yielding salinity, temperature and pressure every 2 m in the vertical from 
the surface and four additional points with salinity and temperature data only at the surface. 
This data was collected between 21 and 25 May off Point Arena, California. Satellite data 
was obtained from channel four of the AVHRR for day 22 May. 

The hydrographic data was used in a first analysis to verify the information contained 
in the satellite image (Figure 4) using three dimensional visualization (section IU.B.l). This 
image shows an eddy of cold water with some warm water possibly entrained from its outer 
edge. The hydrographic data verified the existence of the eddy, and revealed its 
manifestation in the vertical, but showed no clear evidence of the origin of the warm water 
entrainment, possibly due to the low resolution of the data. 

The satellite data was then processed for future use. Since this data is affected by 
atmospheric effects and by a time bias relative to the CTD stations, an adjustment was 
performed to reduce the difference between the satellite temperatures and the CTD surface 
temperatures in a optimal way (section EII.A). 

The hydrographic data was then interpolated from the observations domain into a grid 
domain that can eventually be used in a coastal circulation model (section III.B.2). To 
interpolate the data the method of multiquadric interpolation was used. This method was 
applied in both two dimensional and three dimensional versions, to compare the surface 
results one against the other and against the satellite image. Several runs to make a rough 
estimate of the method’s optimal parameters in the three dimensional case were performed 
(section IV.B). The results obtained are described in chapter IV, with pertinent figures 
contained in the appendix. 
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Based on the research outlined above, the following conclusions are made: 

• The hydrographic data verify the satellite observations, demonstrating that the 
satellite image indeed depicts an oceanic eddy, and not, for example, an 
atmospheric eddy. A particular feature of the eddy, the warm pocket of water 
visible in the satellite image, was also sampled by the CTD casts; 

• This eddy has a vertical signature down to approximately z=-300 m; 

• During the period covered by the cruise the eddy did not undergo any visible 
translation movement thus making the lack of synopticity between the satellite- 
derived and in situ hydrographic data not significant for the purpose of 
objective analysis; 

• The method of muluquadric interpolation does a very good job in extending the 
observed temperature field outside the spatial domain defined by the 
hydrographic data, both in the horizontal and in the vertical; 

• The sensitivity of the multiquadric interpolation to the interpolation constants 
and to the resolution of the observed data is an important factor, particularly in 
the smoothing of the interpolated field. 



During this research some aspects of the data interpolation method suggest work to be 
done in the future: 

• Parameter sensitivity study. This will allow a more accurate determination of 
the optimal parameters; 

• Reduction of small scale noise in the deeper levels of the interpolated 
temperature field; 

• Static stability analysis. The combination of temperature and salinity 
observations to compute the stability of the entire grid domain, thus providing 
a better reference for the performance of the multiquadric interpolation; 

• Inclusion of satellite-derived surface temperatures in the observations vector, 
to improve the description of the temperature field not only in the horizontal 
but also, by means of the multiquadric analysis, in the underlying column of 
water. 
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APPENDIX - FIGURES 
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Figure 1 - Position and number of CTD stations. 
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Figure 2 - Analysis of Sea Surface Temperature (SST). Top figure: comparison of CTD 
and satellite derived SST (SAT); middle figure: difference between the two SST series; 
bottom figure: comparison of CTD and adjusted satellite derived SST. 
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Figure 3 - National Oceanic and Atmospheric Administration (NOAA) satellite image, 
obtained from Channel 4 of the Advanced Very High Resolution Radiometer (AVHRR), 
of the Sea Surface Temperature (SST) off Point Arena, on 21 May 1993, 23 GMT. The 
crosses represent the positions of the CTD stations (temperatures in degrees Celsius). 
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Figure 4 - National Oceanic and Atmospheric Administration (NOAA) satellite image, 
obtained from Channel 4 of the Advanced Very High Resolution Radiometer (AVHRR), 
of the Sea Surface Temperature (SST) off Point Arena, on 22 May 1993, 23 GMT The 
crosses represent the positions of the CTD stations (temperatures in degrees Celsius). 



29 



Latitude 




Longitude 



Figure 5 - Acoustic Doppler Current Profiler (ADCP) instantaneous horizontal velocity 
vectors for depth bin 15-23 meters, for the period 19 to 26 May 1993. 
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Figure 6 - Three dimensional visualization of raw CTD data on a 5x5 station grid. The fig- 
ure shows an isothermal surface (8 °C), several isothermal contours and a vertical section 
of temperature along CTD stations number 5, 6, 15, 16 and 25 (temperatures in degrees 
Celsius). 
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Figure 7 - Three dimensional visualization of raw CTD data on a 5x5 station grid. The fig- 
ure shows an isothermal surface (8.5 °C), several isothermal contours and a vertical sec- 
tion of temperature along CTD stations number 15, 14, 13, 12 and 11 (temperatures in 
degrees Celsius). 
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Figure 8 - Two dimensional visualization of Sea Surface Temperatures (SST) on a 5x5 sta- 
tion grid. The figure was obtained by averaging the first 5 m of raw CTD data (tempera- 
tures in degrees Celsius). The lower left comer corresponds to station number 25. 
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Figure 9 - Grid layout used in the multiquadric interpolation, in UTM coordinates. The X-marks are control points that define 
the grid domain. The crosses represent the positions of the CTD stations. The circles represent the positions of the SAIL sta- 
tions. The origin of the grid (point 0-0) corresponds to geographical position 38.7 ° N, 124.9 ° W. The grid is dimensionless. 
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Figure 10 - Map of SST obtained with the two dimensional version of multiquadric interpolation. The lower left comer corre- 
sponds to the origin of the grid. Red coming from green is cold water, red coming from blue is warm water. 




Figure 11 - Map of Sea Surface Temperatures (SST) obtained with the three dimensional 
version of multiquadric interpolation. The lower left comer of the image corresponds to 
the origin of the grid. Temperatures in degrees Celsius. 
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Figure 12 - Three dimensional visualization of interpolated temperatures, obtained from 
the three dimensional version of multiquadric interpolation. The figure shows a surface 
temperature map, several isothermal contours and an isothermal surface (8 °C). Tempera- 
tures in degrees Celsius. 
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Figure 13 - Three dimensional visualization of interpolated temperatures, obtained from 
the three dimensional version of multiquadric interpolation. The figure shows a surface 
temperature map, several isothermal contours and an isothermal surface (3 °C). Tempera- 
tures in degrees Celsius. 
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